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Abstract 

In the present work the dynamics of a continuous inextensible chain is studied. The chain is 
regarded as a system of small particles subjected to constraints on their reciprocal distances. It 
is proposed a treatment of systems of this kind based on a set Langevin equations in which the 
noise is characterized by a non-gaussian probability distribution. The method is explained in the 
case of a freely hinged chain. In particular, the generating functional of the correlation functions 
of the relevant degrees of freedom which describe the conformations of this chain is derived. It 
is shown that in the continuous limit this generating functional coincides with a model of an 
inextensible chain previously discussed by one of the authors of this work. Next, the approach 
developed here is applied to a inextensible chain, called the freely jointed bar chain, in which 
the basic units are small extended objects. The generating functional of the freely jointed bar 
chain is constructed. It is shown that it differs profoundly from that of the freely hinged chain. 
Despite the differences, it is verified that in the continuous limit both generating functionals 
coincide as it is expected. 

1 Introduction 

There are several physical situations in which it is required that a polymer chain is inextensible. This 
happens for instance when the chain is pulled at its ends by strong forces [H El El HI [5] in experiments 
of DNA micromanipulations. In order to impose the constraints in the stochastic equations governing 
the motion of the chain it is possible to apply a variety of powerful techniques |6l [71 El El [ini ttH [I2l [I3] . 
Progresses in the understanding of statistic and dynamical aspects of inextensible chains have been 
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obtained mainly with the help of numerical simulations, see for example Refs. [HI [151 ttSl [13 UHl [13 
[201 [21] • For analytical calculations, models in which the condition of inextensibility is relaxed are 
usually preferred. This is the case of the gaussian chain [22], which often replaces the freely jointed 
chain [23] in the study of the statistical behavior of polymers systems. Other examples of theories in 
which the length of the chain is allowed to become infinitely long are the models of Rouse [21] and 
Zimm [25]. They provide a satisfactory description of the dynamics of polymer chains in solutions, 
but could fail in the case of a chain that is stretched under strong external forces applied at its ends. 

For these reasons it is highly desirable to construct a path integral model describing the dynamics 
of an inextensible chain which is simple enough to allow analytic calculations. In particular, we are 
interested in the derivation of the generating functional of the correlation functions of the relevant 
variables that describe the conformation of the chain. The usual starting point is a discrete mechanical 
model which provides a coarse grained approximation of the chain. Successively, a continuous limit 
is performed, in which the number of elements becomes infinitely large and their length approaches 
zero, while the total length of the chain L is kept constant. This limit is useful in order to smooth 
out the dependence on the details of the particular mechanical model chosen. In the first part of this 
article the case of a freely hinged chain (FHC) is considered. This is a system of beads connected 
together by massless links of fixed length a. The FHC has been studied in the past in Refs. [2S1 [23 [2S] 
in connection with the dynamics of cold polymers or a gas of hot polymers. In Ref. [21] it has been 
shown that the generating functional for a FHC coincides in the continuous limit with the partition 
function of a theory which closely resembles a nonlinear sigma model. For this reason it has been 
called the generalized nonlinear sigma model or simply GNLaM. Within the GNLaM it is possible 
to perform analytical calculations of physical quantities that can be applied to study the dynamics 
of cold chains or of chains moving in a very viscous environment [29l [30l |3T]. One drawback of 
the GNLaM is that it regards the FHC as a gas of fluctuating particles with the inextensibility 
constraints being implemented by means of a product of delta functions. The introduction of delta 
functions to fix constraints is a standard procedure in the investigations of the statistical mechanics 
of polymers, where it is an useful tool in order to limit the conformations of the chain. In the case 
of dynamics, however, a physical explanation of the appearance of the delta functions inside the 
path integral of the GNLcxM is necessary, otherwise the connection with the stochastic process of the 
fluctuating chain is lost. Such connection was up to now somewhat obscure, despite the fact that 
in Ref. [32] it has been possible to verify that the GNLcrM reproduces as expected the sum over all 
chain conformations satisfying a free Langevin equation and subjected to constraints that ensure the 
property of inextensibility. 

One goal of this work is to provide a physical interpretation of the GNLcrM as a stochastic process 
in which the noise is non-gaussian. The starting point is the observation that the chain degrees of 
freedom which are constrained by the inextensibility requirement are frozen and thus cannot be 
influenced by the noise. This implies that the noise must be constrained too. The conditions on 
the noise are imposed in the present approach in a soft way using an elastic potential. It is the 
introduction of this potential that modifies the noise distribution and makes it non-gaussian. The 
rigid constraints are recovered in the limit in which the elastic constants become infinite. Following 
this strategy, the motion of the beads of the chain is described by a system of free Langevin equations 
in which the random forces are characterized by a non-gaussian distribution. For any finite value 
of the elastic constants, the related generating functional corresponds to that of a chain consisting 
of beads connected together by springs. With increasing values of the elastic constants, the springs 
become stiffer and stiffer until at the end an inextensible chain is obtained. Let us note that also 
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the beads of the Rouse model are joined together by springs. However, in that model the potentials 
are chosen in such a way that the springs are in their rest positions only when the distances between 
neighboring beads is zero. For this reason, in the limit of infinite elastic constants, the chain collapses 
to a point. Other differences of the GNLaM from the Rouse model have been discussed in details in 
Ref. [29]. The present strategy to enforce the inextensibility constraints with the help of potentials 
closely resembles the way in which stiff constraints (sometimes called flexible constraints) are imposed 
[201 EI] in stochastic equations. Stiff constraints are fixed in fact with the help of potentials [331 [31] 
whose role is to create an energy barrier that strongly suppresses those chain conformations that do 
not satisfy the constraints. The form of these potentials, like for instance the FENE potential, is 
usually too complicated to allow analytic calculations. In the case of the GNLaM, instead, there are 
only elastic potentials, which are relatively simple. Another striking difference with respect to stiff 
constraints is that in our approach it is taken the limit in which the strength of the elastic forces 
become infinite. In this limit the constraints become rigid and not flexible. 

In the second part of this work the idea of Ref. [iSj is investigated. According to that idea, a 
chain may be realized by connecting together basic inextensible units, like for instance small bars. In 
this way, the length of all elements of the resulting chain is always constant by construction without 
the need of imposing the cumbersome holonomic constraints that are otherwise necessary to enforce 
the property of inextensibility. Constraints are still necessary in order to connect together the basic 
units, but they are very simple and their imposition is straightforward. In deriving the generating 
functional for a system of this kind, however, the fact that the basic inextensible units are not point- 
like objects, but rigid bodies with rotational degrees of freedom, is the origin of several complications. 
For this reason, in this work the basic inextensible units are built starting from a number K of beads 
of diameter A/, called hereafter spheres. The motion of the spheres is constrained in such a way 
that they stay aligned along a segment of fixed length. The shape of the whole system formed 
by the K spheres is that of a shish-kebab and offers a good approximation of a bar. Let us note 
that the constraints needed in order to form the shish-kebabs consist in conditions on the distances 
between pairs of points. Thus we can apply the same strategy used to implement the inextensibility 
constraints in the FHC, which are of the same type. In the limit in which the number of spheres 
approaches infinity and their radius becomes vanishingly small, while the length of the basic units 
remains constant, one obtains what has been called here the freely jointed bar chain or simply FJBC. 
This is a discrete chain composed by one- dimensional bars with uniform distribution of mass. 

The expression of the generating functional of the correlation functions of the shish-kebab chain 
has been provided in the most general case of a chain containing basic units, each composed by a 
number K of beads. The generating functional of the FJBC in the limit K — >■ +infty and A/ — )■ 
has been derived. As discussed in [5] , the constraints that are needed to held together the basic units 
are almost trivial and the associated spurious degrees of freedom can be easily eliminated. Finally, 
the continuous limit of the generating functional of the FJBC has been performed. In this case several 
simplifications occur and, at the end, it has been possible to prove that the generating functional of 
the resulting continuous system coincides with the GNLaM. This is an expected result, because in 
the continuous limit the mechanical details of the underlying discrete chain should not play any role. 

The material presented in this paper is organized as follows. In Section [2] the classical aspects of 
the FHC and the FJBC are analyzed. For the interested reader, a more extensive discussion on this 
subject can be found in Ref. [35]. First, the case of the FHC is briefly reviewed. Next, the FJBC is 
studied. It is shown that, while the kinetic and potential energies of the FJBC differ considerably from 
those of the FHC, they coincide in the continuous limit. The discussion of the classical aspects of the 
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chain dynamics is very important to establish a mathematically convenient description of the FJBC, 
that will be used later to study its statistical dynamics. The classical dynamics of inextensible chains 
is also very interesting for concrete applications ranging from cosmetics to computer graphics. Section 
[3] is dedicated to the statistical dynamics of the chain. In the first part of that Section, the generating 
functional of the correlation functions of the radius vectors which specify the positions of the beads 
of the FHC is derived in path integral form. The constraints required by the fact that the length 
of the links connecting the beads is constant are imposed by using elastic potentials as explained 
before. When the strength of these potentials becomes infinite, the path integral formulation of the 
generating functional of the FHC given in Ref. [29] is recovered. In the continuous limit the GNLcxM 
is obtained. We fill also a gap of our previous publications concerning the FHC by establishing the 
connection of the partition function of the FHC with a stochastic process described by a Fokker- 
Planck equation. In the second part of Section [3] the generating functional of the FJBC in the 
shish-kebab approximation is derived. Next, in Section H] we compute the generating functional of 
the FJBC without any approximation. In analogy with what happens in the classical case, it turns 
out that this generating functional is very different from that of the FHC. However, we are able to 
show that they both coincide in the limit of a continuous chain. Finally, our conclusions are drawn 
in Section O 



We begin by studying the dynamics of a chain composed by beads of mass m joined together by 
— 1 massless segments of length a. Though it is not strictly necessary, to avoid complications with 
the spherical coordinates in arbitrary dimensions, in this Section the discussion will be limited to two 
dimensional chains moving on the {x,y) plane. Let Ri(t) = {xi(t),yi(t)) denote the radius vector of 
the i—th bead, for i = 1, . . . , A^. In polar coordinates (see Fig. [1]): 



2 Classical dynamics 



2.1 The freely hinged chain (FHC) 



i-l 




I = 



(1) 




I = 



2,... 



N 



(2) 



We suppose that the first point {i 



1) is fixed in the origin, so that: 



x,{t) = yi{t) = 



(3) 



First of all, we compute the kinetic energy of the system: 



N 




(4) 



i=2 



4 




y 



Figure 1: FHC consisting of N beads of mass m (the black circles) connected together by — 1 
segments of fixed length a. The angles between the y axis and the i—th segment are denoted by a^, 
i = l,...,N -1. 

The dots denote here derivatives with respect to the time t, i. e. x = After a few calculations we 
find the expression of Tfhc{o) in polar coordinates: 

2 i-1 

Tfhc{o) = ^ ^ ajcxk cos(aj - a^) (5) 

i=2 j,k=l 

We are interested in the continuous chain obtained from the FHC after performing the limit in which 
the length a of the segments goes to zero while the total length L of the chain is preserved: 

a — ^ — y +00 and Na = L (6) 

In that limit the sums become integrals according to the prescription: 




A function fi{t) depending on the discrete index i, like for instance the angles ai{t), will be substituted 
by a function of the continuous variable s, where < s < L, according to the rule: fi{t) — > f{t, s). 
The full procedure to pass to the continuous limit is described in Ref. [23 ES] and will not be 
discussed here. The new parameter s appearing in ([7]) is the arc-length that measures the distance 
between two points on the chain. Partial derivatives with respect to s will be denoted with a prime, 
for instance R'{t, s) = ^^j^- We also need to express the mass of the beads m as a function of a. To 
this purpose, we introduce the density of mass p on the chain. Since the mass distribution is uniform 
- all beads have the same mass - we may put P = j;, where M is the total mass of the chain. Thus, 
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m can be written function of the mass density: 

M 

m = —a 
ij 



It is easy to check that, for consistency, Nm = M because Na = L. Applying the above prescriptions, 
it is 

by: 



it is possible to show that the continuous kinetic energy TpHc = 1™ TpHcio^) is given 

a^O,N^+oo,Na=L 



7f//c = ^ / ds I du I dva{t,u)a{t,v) cos {a{t,u) — a{t,v)) (9) 







Let us note that the expression of TpHC given above is different from that of Ref . p9] . This is due to 
the fact that here we have not followed the recursive procedure for computing TpHc{(^) used in [29] . 
but rather we have derived it directly from the definition of kinetic energy of Eq. (jl]). It is possible 
to show the equivalence of the results obtained in this work and in Ref. |29] exploiting the identity: 



du I dvf{u,v)= du dvf{u,v)+ du dvf{v,u) (10) 



which is valid for an arbitrary integrable function f of v and u. Exploiting Eq. ( ITOj) . the kinetic 
energy TpHc in © becomes: 



M 



TpHC = -j- / ds du dva{t,u)a{t,v) cos{a{t,u) — a{t,v)) (11) 
L Jo Jo Jo 

The right hand side of the above equation coincides exactly with the kinetic energy of the FHC 
derived in |29j . 

Finally, we add to the FHC model the interactions. To this purpose, we define the potential: 

AT N 

VpHcia) = Y,(^VE{R,{t)) + a^Vj{R,,{t),Rj{t)) (12) 
1=1 «j=i 

where Vp^Ri^t)) and Vi{Ri(t), Rj{t)) have the meaning of energy densities per unit of chain length 
due to external and internal forces respectively. For this reason, in the total potential energy of the 
chain VpHc{(^) appearing in Eq. (IT^ . Vp{Ri{t)) and Vj{Ri(t), Rj(t)) are multiplied by the factors 
a and respectively. The presence of these factors is in agreement with the continuous limit 
prescription of Eq. ([7j). 

In the continuous limit ([6]) we obtain from VpHcio-)'- 

VpHc= lim VpHc{ci)= dsVp{R{t,s))+ ds ds'Vi{R{t, s), R{t, s')) (13) 

a^O,N^+oc,Na=L Jq 

where R{t, s) is the continuous limit of the bond vectors Ri{t). From Eqs. and ([2]) it turns out 
that the components of R{t, s) are given by: 

R{t,s)= / du{sma{t,u), cos a{t,u)) (14) 
Jo 
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2.2 The case of the FJBC 

Let us consider a chain composed by N bars of length a. Each bar is regarded here as a one- 
dimensional segment with a uniform distribution of mass along it. The bars are joined at the points: 



i-l 
' =1 



= ^^asina-' . (15) 



i-l 

yi = ^ a cos a/ . (16) 

i'=i 

where i = 2, . . . , N . One end of the first bar is supposed to be fixed in the origin, so that Xi — yi — 0. 
The coordinates of every point of such a chain are given by: 

Xi {t, I) = asin aj{t) + / sin ai{t) (17) 

i-l 

Hi {t, I) = a cos aj{t) + / cosai(t) (18) 

ior i — 2, . . . , N — 1. The points of the first bar have coordinates: 

Xi {t, I) — I sin ai{t) (19) 

yi{t,l)^ I cos ai{t) (20) 

Here I is a continuous parameter taking its values in the range < Z < a. In the future also the 
notation 

Ri{t,l)^{xi{t,l),yi{t,l)), i^l,...,N-l (21) 
will be used. We are now ready to compute the kinetic energy of the FJBC: 

T^jBcia) = £ ^dl {xUt, I) + yKt, 0) (22) 

In the above equation /li is the mass of the i—th bar. We have assumed that the mass distribution 
of each bar is homogeneous. For that reason, the mass dnii of a small element of length dl is: 

drrii = —dl (23) 
a 

If all bars have the same mass, i. e. /Xj = /x for i = 1, . . . , — 1, it is possible to put drrii = ^dl 
using the fact that /iN — M and aN — L. After simple calculations we find out that: 



■ 1 «/ 



i-l i-l 

a^ajCXk cos(Q;j — a^) + Pa'^ + 2al djCtj cos(Q;j — ctj) 

1 \j,k=l 3=1 



(24) 
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After the integration over dl we obtain: 



3 

i=l 



i—1 . 2 1 



(25) 



ajCtk cos(aj — a^) + + ctjOa cos{aj — ai 

.j,k=l j=l 

Let us note that, in order to derive the above expression of the kinetic energy, the computation of the 
moments of inertia of the bars has not been necessary. Each bar has been rather decomposed into 
a set of infinitesimal elements which can be regarded as points. The kinetic energy of each of these 
points has then been summed up to obtain the total kinetic energy of the chain given in Eq. ( l25l) . The 
result is the same as if we had considered the bars as the basic units of the chain and computed their 
kinetic energy. An approach in which the chain is regarded as a set of points and not as a system of 
bars is not a great advantage in the simple classical case. However, it will be of crucial importance 
when the statistical physics of the FJBC will be studied, because the Fokker-Planck equation for a 
bar is prohibitively complicated for our purposes. 

We are now ready to perform the continuous limit (E]) in which the length a of the segments goes 
to zero while the total length L of the chain is preserved. Putting 

lim ^ TpjEcia) = TpjBC (26) 

a^O,N^+oo,Na=L 

the continuous kinetic energy takes the form: 

M r r 

TpjBC = —f I ds du dva{t,u)a{t,v) cos{a(t,u) — a(t,v)) (27) 
2iv Jo Jo Jo 

Eq. (1271) coincides exactly with the kinetic energy of the continuous FHC of Eq. ([9]) . 

Let's now add to the FJBC the interactions described by a general potential of the kind: 

VFjBc{a)=J2 dlVE{R^{t^))+Y, dl dl'Vi{R,{t,l),R,{t,l')) (28) 

*=1 i,j=l -^0 ^0 

Here VE{Ri) and Vi{Ri, Rj) take into account the external and internal interactions respectively. In 
the case of small values of a, Vfjbc («) may be approximated as follows: 

N-1 N-1 

VrjEcia) ~ ('^EiRiit, a))+ ^ a^Vi{Ri{t, a),R,{t, a)) (29) 

i=l i)i=l 

with 

i-l 

Ri{t,a) ='Y^a{ sin a j(t), COS a j{t)) i = 2,...,N (30) 

i=i 

and Ri{t) = 0. The above approximation of Vfjbc'{0') becomes exact in the continuous limit of 
Eq. (E]), in which the expression of the potential becomes: 

pL pL pL 

Vfjbc = dsV E{R{t, s))+ ds ds'Vi{R{t, s), R{t, s')) (31) 
Jo Jo Jo 

where R{t,s) is given in Eq. (fl^ . As we may see, the above expression of the potential coincides 

with that of the FHC of Eq. (fT3|) . As a consequence, in the continuous limit the FHC and FJBC are 

equivalent as expected, despite the fact that the discrete models have different kinetic and potential 

energies. 
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3 Chain statistical dynamics 



3.1 The case of the FHC 

In this section we consider the motion of a FHC composed by beads in d dimensions. Each bead 
is fluctuacting in a viscous medium kept of constant temperature T. The viscosity makes the motion 
overdamped [37J, so that the radius vectors of the beads satisfy the Langevin equation: 

Ri = Ui{t), t = l,...,N. (32) 

where Vi{t) is a random force with a probability distribution that will be specified later. 
Due to the constraints, the above equations have to be completed by conditions: 

a{Ri, R^-i) = \Ri{t) - Ri-i{t)\ - a = 0, t = 2,...,N. (33) 

It is easy to solve Eq. fl32l) with respect to the Ri's. The result is: 

Ri,uAt)= [ iyi{r)dr + Ro^i (34) 

We have denoted the solutions of Eq. (1321) with the symbol Ri^i,^{t). The superscript h>i is used to 
stress the dependence of these solutions on the noise. Of course, not all components of the noise 
are independent. In fact, if we plug in the solutions -Rj,y.(t)'s of Eq. fl32|) in Eqs. fl33l) . we obtain 
constraints for the noises i^j's: 

a(i?i,,.,, = \Ri,uAt) - R^-l,u.-,it)\ - a = 0, i = 2,...,N. (35) 

The system of equations fl35|) allows in principle to eliminate the degrees of freedom of the random 
forces i>i that become spurious because the beads are fixed at the ends of segments of constant length 
a. Unfortunately, even in the present simple case in which the external forces are absent, it is not 
easy to solve Eqs. (l35l) by expressing the redundant — 1 degrees of freedom as a function of the 
remaining 2N + 1 independent variables. 

An alternative procedure consists in the introduction of a set of Lagrange multipliers Aj(t), i = 
2, . . . , N in order to impose the constraints. In this approach the i/j's are regarded as unconstrained 
sources of white noise. The price to be paid for that is the addition of the reaction forces 

dCi{Ri,Ri^i) . 
FR,i = Xi — i = 2,...,N (36) 



dR^ 



in the Langevin equations (1321) . The solutions Ri^i^ x. of the new Langevin equations obtained in this 
way will depend also on the Lagrange multipliers A,. The latter can be determined by exploiting the 
constraints ( 133|) . In practice, it is more convenient to use the consistency relations coming from the 
requirement that the constraints should be preserved in time: 

Ci(-Ri,i.,,A,, J?i-i,iy,_i,A,_i) = (37) 

Explicitly, we have that 

C'j(-Rw,Ai, -Ri-i,i/,_iAj = -{,Ri,u„\, - Ri-i,u,.^,\,) ■ {i^i - J^i-i) = (38) 
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for i = 2, . . . , N . The constraints fixed with the help of the Lagrange multiphers are sometimes called 
rigid constraints. The resulting Langevin equations are mathematically complicated, but yet their 
solutions may be derived numerically. 

With both methods, direct solution of Eqs. f p5|) or the introduction of the Lagrange multipliers, 
it is very difficult to arrive to a path integral formulation of an inextensible chain that can be used 
to perform analytical calculations. We recall that we are interested here in the construction of the 
so-called generating functional of the correlation functions of the solutions of the Langevin equations 
-Rj It will be shown in this Section that it is possible to achieve this goal starting from a slightly 
different point of view. The idea is to replace the usual gaussian noise distribution by the probability 
distribution Vp{{vi}) given by: 



N 



exp 



N 



AD 



(39) 



The potential V^(i?i^^^, . . . , Rn^^^, (32, ... , [3n) will be selected according to the following criteria: 

1. When P2 = Pi = ■ ■ ■ = Pn = ^-1 the system is unconstrained and the i/j's become purely 
gaussian noises. 



2. For large values of the /3i's, the potential V{Ri 



R 



, Pn) should exhibit a sharp 



minimum near the region of noise configurations in which the constraints (135!) are satisfied. 



3. Finally, in the limit /?, = +00, the potential V{Ri^i,^, 
if Eqs. (135 p are not fulfilled and zero otherwise. 



, RN,iyN, P2, 



(32, ■ ■ ■ , Pn) should be infinite 



Clearly, thanks to the third condition all noise configurations that do not conform with Eqs. (!35|) are 
eliminated. Only those for which the constraints ( 135|) are satisfied remain. 

To derive a potential with the characteristics described above we can use the physical intuition 
according to which a chain made of beads and springs will behave as a chain of beads and links of 
fixed length in the limit of infinitely large spring elastic constants. Thus, we consider here a system 
of beads connected together by iV — 1 springs. Let (3i, i = 2, . . . , N he the elastic constant of the 
i—th spring and a its rest length. The elastic forces: 



(3ii\Ri-R,^i\-a) 



[Ri — R, 



\Ri ~ Ri-l\ 

are acting on the internal beads with indexes i = 2 

F,=(32{\Ri-R2 



(3i {\Ri — Ri+i\ — a 



[Ri — R 



■i+l, 



\Ri 



Ri+l\ 

N — 1. The following forces 
(Ri - R2) 



.N 



1) 

(40) 
(41) 

\R R I ^^2) 
\JrCN — rCN~i\ 

are applied instead on the beads lying at the ends of the chain. While from the point of view of the 
mathematical complexity there is no problem in considering — 1 independent elastic constants 
in practice this is an unnecessary complication. For this reason, from now on we will assume that 
the /3i's are all equal: 

(3i = (3 1 = 2,. ..,N (43) 



(3n{\R 



■N 



R 



\Ri- 
(Rn 



R2\ 

- Rn-i) 



N-l\ 
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After this simplification, the potential corresponding to the interactions (H0] - H2|) may be written as 
follows: 



N 



V{Ri, ...,RN,P) = ^aVj{\Rj- R,.i\, 
i=2 



/3) 



where 



V,{\R,-R,^,l^) = ^i\R, 



R 



(44) 



(45) 



The minima of the potentials (H5l) and thus the minimum of V{Ri, . . . , Rn, P) occur when the 



conditions \Rj — Rj-i\ — = 0, j 



N, are verified. These are exactly the inextensible 



constraints of Eqs. f l33|) . It is easy to show that V{Ri, . . . , Rn,P) satisfies requirements 1. - 3. 
We are now ready to write down the expression of the generating functional of the FHC, which 



is defined as the average of the quantity exp — a /o^ '^j(^) ■ Ri,ui{t)dt 
probability distribution Vp{{vi}) of Eq. ( 139|) : 



with respect to the 



'FHC 



[J] 



lim 



W vum exp -/ * E^+E^fil 

,1=1'' J [ "'0 \i=l j=2 



R 



R 



,/3 



exp 



N 
i=l 



J,{t) ■ Ri,^Xt)dt 



(46) 



Let us note that in Eq. (146|) the path integration is extended over all possible noise configurations. As 
already anticipated before, the constraints are fixed using the potentials Vj (|-Rj>j — Rj-i,uj^i \ , 1^) 
after taking the limit /3 — > +oo. Indeed, when /3 is very large all configurations Ri^ui that do 
not satisfy the conditions Ci{Ri^^-, Ri_i^^-_J = are exponentially suppressed. Eq. (H^ should be 
completed by specifying the positions Ri^u, of the beads at the initial and final instants. We require 
that i2j j^^(0) = i^o.j according to Eq. (^^. Moreover, when t = tf it is assumed that the z— th bead 
is located at the point Of course, both iio,j and Rf^i should satisfy the constraints fl33|) . 
At this point we insert in Eq. ( H6ll the quantity 



N 

n 



(47) 



^ Jr,{0)=Ro,, 



The boundary conditions of the integration over the new fields T*i(t) have been chosen in such a way 
that they are consistent with the boundary conditions of the fields Ri^u,- Clearly J = 1, so that the 
insertion of / in Eq. (146|) will not change the physics of the problem. As a consequence, remembering 



the explicit expressions of the potentials Vi{\Ri 
follows: 



R 



i—l 1 5 



/3) of Eq. ( H5|l . we may rewrite ZpHc[J] as 



'FHC 



[J] 



lim 



■ TV 

n 



2 = 1 



exp 



N 



i=l 



J,{t) ■ n{t)dt 



exp 



exp 



N 



i=l 



N 



,i=l 



3=2 



of dt 



(48) 
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Now we recall the identity: 



6{ri - Ri^^J = det 



d_ 

di 



6{ri - Ui 



(49) 



N 

where we have used the fact that Ri^^i satisfies Eq. fl32|) . Ignoring the constant factor Y[ det 

1=1 



and performing the simple integrations over the noises i/,, the generating functional Zfhc[J] takes 
the form: 



Zfhc[J] 



lim 

/3-J-+00 



X exp 



■ N 

n 



ri{tf)=Rf^i 



^ Jri(0)=Ro,, 



exp 



i=l 



3=2 



(50) 



Before continuing, a digression is in order. First, we rewrite Zfhc[J] as follows: 

Zfhc[J] = lim [vpp{{r,{t)})e-^^-^^o'-J^-r. 
where Vpp{{ri(t)}) is the probability distribution: 



(51) 



N 



i=l 
N 



exp 



X 



n[5(r-.(t/)-i2/,)5(r-.(0)-ilo,i)] 



(52) 



i=l 



The last two delta functions are needed to impose the boundary conditions. For each fixed value of (3 
the probability function Pg({_Rj^j}, tj; {i?o,i},0) of the system of beads and springs is defined as the 
integral over all possible configurations of the above probability distribution [38] : 



Pp{{Rf,},tf;{Ro,},0) 



■ N 

n 



r,{tf)=Rf^, 



1 Jr,{0)=Ro,, 







exp 









N 



1=1 



The probability function P{{Rf^i},tf; {Ro^i},0) of the FHC is obtained in the limit (3 — )■ +oo: 

Pi{Rf,}, tf, {RoA. 0) = lim Pp{{Rf,i}. tf, {RoA, 0) 



(53) 



(54) 



Let us note that P^({i?/ i}, t/; {-Rq,*}, 0) satisfies the Schrodinger-like equation of a system of particles 
with interactions described by the elastic potential ( l44l) : 



N 



■I 1 = 1 



P^({i?/,},t/;{i?o,},0) = 



(55) 



As it is possible to see, the elastic interactions do not appear inside a drift term as it happens in 
standard Fokker-Planck equations. This is not a surprise, because the elastic forces do not describe 
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any physical property of the chain. They have been introduced in the non-gaussian noise distribution 
of Eq. (15^ with the sole purpose of imposing the constraints. On the other side, the generating 
functional Zfhc[J] is related to a stochastic equation. To show that this is exactly the case, we use 
the well known connection between Schrodinger-like and Fokker-Planck equations. First of all, we 
introduce the new potential U (-R/,i, • • • , Rf,N, P) through the differential equation: 



D 



N 

E 

i=l 



9Rh 



[Rf,i, • • • 5 Rf,N, (3) + 



dU 



OR 



■Li 



{Rf,i, . . . , Rf,N, (3) 



V{Rf,,,...,Rf,N,/3) (56) 



Let us note that U{Rf^i, . . . , Rj^n, P) is dimensionless, so that we may rescale the partition function 
Pg({i?/,i},t/; {i?o,i},0) as follows: 



P^({i?/,,},t/;{i?o,},0) = P^({i?y,},ty;{i?o,.},0)e^(^^'i'---'^/'-'^^ 



(57) 



It is easy to check that the new partition function P'p{{Rf^i},tf] {-Ro,i}) 0) satisfies the Fokker-Planck 
equation: 



_d_ 



1 ^ 
5E 



d 



where the forces fi{Rf^i, . . . , Rf^N, P) are given by: 

fi{Rf,i-, ■ ■ ■ , Rf,N, /3) = —2D 



d 

^TTB ^ fi{Rf,l^ ■ ■ ■ ; Rf,N, /3) 



P' {{Rf,},tf,{Ro,}, 0) = (58) 



^UiRf,^,...,Rf,,,,/3) 



OR 



■Li 



(59) 



This concludes our proof. 

Let's now go back to the generating functional of Eq. (l50l) . Its expression may be further simplified 
exploiting the relation: 



(60) 



lim e-^'^^o^ f{r,{t),r,{t))dt ^ S{f{ri{t), rj{t))) 

which is valid up to an (infinite) proportionality constant for a generic function f{ri,rj). Eq. ( l60l) 
will be proved in Appendix A. Applying Eq. ( I60l) in Eq. ( 150|) in the limit of large values of (3, we 
obtain: 



'FHC 



[J] 



■ N 

n 



Vri{t) 



3 rj-i\-a) 



exp 



i=l 



' aJj ■ Vj ) dt 



■ N 

.i=2 



(61) 



The above generating functional coincides, apart from an irrelevant constant, to the generating func- 
tional of the discrete chain composed by beads and links discussed in Ref. |29j. To show that, first we 
note that inside Eq. ( 16T|) it is possible to replace the functional Dirac delta functions 5 {\rj — rj_i| — a) 
with 5 {\rj — Tj_ip — a^). A rigorous proof of this fact has been provided in Ref. [29]. Basically the 
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proof consists in the generalization to functional delta functions of the following delta function rela- 
tion 

S{x^ -a^) = ^ {S{x -a) + 6{x + a)) (62) 

and on the consideration that the set of chain conformations for which \rj — rj_i| + a = is empty. 
Using also the identity: 



■ N 

J=2 



(63) 



and neglecting irrelevant constants we may rewrite Eq. ( l6T|) as follows: 



'FHC 



[J] 



■ N 

n 



ri{tf)=Rfj 



^ Jri(0)=i?o,i 



exp 



N 
i=l 



—rj + aJi -Vi] dt 



■ N 

J=2 



(64) 



This is exactly the generating functional of the FHC derived in Ref. [29]. The continuous limit of 
Zfhc[J] is known and coincides with the partition function of the GNLuM: 



Z[J] 



r{tf,s)=Rf{s) 



r{0,s)=Ro{s) 



Vr{t, s)e'^^^^°^ dtf^^dsf^it,s)^-f^f dtf^'^dsj(t,syr(t,s)^ 



dr{t, s) 



ds 



(65) 



Here ks is the Boltzmann constant, T is the temperature and r is the relaxation time of the beads 
composing the chain. These quantities arise from the factor ^ appearing in Eq. and can be 
derived following the procedure explained in [29] and repeated here for convenience. First, we note 
that D = fikT, where /i is the mobility of the beads. Remembering the fact that = ^, we get 
D = and thus ^ = Exploiting the fact that m = ^a, we obtain the desired expression: 



1 



1 M 



AD 2kTT 2L 



(66) 



Z[J] describes the dynamics of a continuous chain of length L. The inextensibility constraints are 
imposed by the functional delta function 6 ( ^^^j^ 



3.2 The case of the FJBC 

With respect to Subsection 12.21 which was dedicated to the classical FJBC, in this Subsection we 
slightly change the definition of the small basic units composing the FJBC. Instead of one-dimensional 
segments with uniform mass distribution, each bar is replaced by a shish-kebab model consisting of 
K small spheres with diffusion constant D and mass m as shown in Fig. [2J The distance Al of each 
sphere is given by: 

Al = ^ (67) 
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Figure 2: Shish kebab model for the bars composing the chain. 



In the case of dynamics it is possible to set the diameters of the spheres by choosing their diffusion 
constants (or alternatively their mobility, their relaxation times etc.) appropriately. Here we will 
assume that the diffusion constant is given by D = 3^^^, where rj is the viscosity of the medium in 
which the chain fluctuates. This choice corresponds to spheres of diameter Al. It does not reduce 
the generality of our discussion, which remains valid in the case of any other choice. The positions of 
the centers of mass of the k—th sphere belonging to the i—th shish-kebab is described by the radius 
vectors Ri^k, i = I, ■ ■ ■ ,N and k = 1, . . . ,K. Throughout the rest of this Section the words shish- 
kebab and bar will be used interchangeably despite the fact that they are related to different objects. 
In the next Section we will see that one-dimensional bars are recovered in the limit K — )■ +00. 

At this point we are ready to impose the constraints that fix the positions of the spheres in such 
a way that their motion will not destroy the shape of the FJBC First of all, we need to derive the 
set of conditions satisfied by the spheres belonging to the same bar. To this purpose we note that if 
two spheres k and k' belong to the same i—th bar, the distance between their centers of mass must 
be I A; — A;' |. As a consequence, the radius vectors Ri^k and Ri^k' have to satisfy the relations: 

\Ri,k'it)-R,,k{t)\ = ^\k-k'\ k^k' k,k' = l,...,K (68) 

The above set of constraints is sufficient in order to guarantee that the spheres in a bar will remain 
aligned during their fluctuations. 

Next, we have to be sure that the ends of the bars are correctly connected together in order to 
form a chain. This goal is realized by identifying the K—th sphere of the i—th bar with the first 
sphere of the {i + 1)— th bar, see Fig. |3l This requirement implies that the locations of the centers 
of mass of the spheres should be constrained by the conditions: 

Ri^K = Ri+i,i i = l,...,N-l (69) 

Let us note that Ri^i and Rn,k are not affected by the constraints (!69l) because they are at the ends 
of the chain which are free. 

Before continuing our main discussion, let's make a digression concerning the length of the FJBC 
defined above. The total length of the chain is not L = Na, because we have to subtract from this 
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Figure 3: This figure shows how the bars are connected together to build the FJBC 



value the length A/(A^ — 1) of the spheres which are identified at the joints. After doing that, we 
obtain the effective length L of the chain: 

L = Na- Al{N - 1) (70) 

Due to the fact that AZ = we get: 

Let us now come back to the constraints (l68l) and fl69l) . They represent all the constraints that 
are needed to make sure that the shape of the FJBC is preserved during the motion of the spheres. 
Luckily, the constraints ( l69l) may be easily eliminated. This fact was already noted in Ref. [5] 
and is one of the main advantages to construct an inextensible chains starting from a set of rigid 
and inextensible basic units. In order to get rid of the conditions ( 169!) . we choose as independent 
coordinates the radius vectors -Rj+i_i's. Accordingly, the radius vectors Ri^K^s have to be replaced 
everywhere with the -Rj+i^i's. First of all, we have to perform the substitutions Ri^x — > Ri+i,i 
inside the relations fl68p . In order to proceed, it is convenient to divide these relations into three sets. 
The first set imposes conditions only on the coordinates Ri^k for i = 1, . . . , N and k = 1, . . . , K — 1. 
These coordinates are not located at the junctions between two neighboring bars and thus are not 
affected by the elimination of the constraints fl69|) . As a consequence, in the restricted range of the 
indexes for which i = 1, . . . ,N and k ^ k' = 1, . . . , K — 1, we can simply rewrite the constraints fl68l) 
without changes: 

\Ri,k' - Ri,k\ = T^\k - k'\ k^k' kj^k' = l,...,K-l z = l,...,N (72) 

Let's consider now Eq. (|68l) for k' = K , k = 1, . . . , K — 1 and i = 1, . . . , N — 1. This set of conditions 
contains the coordinates Ri^x for i = 1, . . . , N — 1. They should be substituted by the variables 
R i+i^i's as mentioned above. This operation results in the new conditions: 

\R^+l,l - R^,k\ = -^\K - k\ t = l,...,N-l k = l,...,K-l (73) 
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At this point, we are left only with the subset of equations ( 168|) which is strictly related to the spheres 
of the A^— th bar. Clearly, the K—th sphere in the iV— th bar is on one of the two free ends of the 
chain. Thus, it is not constrained by Eq. f l69l) and the constraints of Eq. (1^5]) should not be changed 
in this case: 



\R 



N,K 



R 



N,k\ 



a 
K 



K-k\ 



k 



K 



(74) 



We are now ready to proceed with the construction of the generating functional of the FJBC Ex- 
actly as we did for the FHC, we first introduce the Langevin equations that describe the fluctuations 
of the spheres. In the present situation the number of spheres is KN . Their positions are denoted 
by the vectors Ri^ki where i = 1, . . . ,N and k = 1, . . . , K. Accordingly, we have to introduce a set of 
KN Langevin equations: 

Ri,kit) = lyi^t) t = l,...,N k = l,...,K (75) 

In analogy with the case of the FHC, the probability distributions of the noises Vi^^s are nonlinear 
as a consequence of the constraints of Eqs. (172147^ . Before giving the explicit expressions of these 
distributions, however, we have to impose in equations fl75l) the constraints fl69l) . Due to these 
constraints, in fact, not all the random forces appearing in Eq. fl75|) are independent. Indeed, using 
the relations fl^ it is easy to check that: 

i = l A^- 1 



(76) 

The above equations have a very simple physical explanation. In order to join the bars together, 
some of the spheres have been identified. As a result, the spheres identified in this way are also 
subjected to the same noise and thus Eqs. (!76|) should be satisfied. In a similar way as we did for 
the constraints (1681) . we may easily get rid of the redundant noises fi,Ki i = 1, . . . , N — 1 from the 
Langevin equations (!75l) at the price of dividing these equations into two sets. In the first set there 
will be only the spheres with indexes k < K: 



N 



k = l,...,K -1 (77) 
■ 1 are spurious due to the constraints 



Ri,k — l^i,k ^ — 1, 

The degrees of freedom Ri^K and Ui^x for i = 1, . . . , N 
and (I76|) . so that they should not be taken into account. It remains only to consider the fluctuations 
of the last sphere on the last bar corresponding to the indexes i = N and k = K. This sphere is not 
constrained because it is located at one free end of the chain. The related Langevin equation reads 
as follows: 

Rn,k = l^N,K (78) 
There are no other independent Langevin equations besides Eqs. fl77|) and f l75]) . 

Having eliminated the constraints ( 169|) . we are ready to construct the generating functional of the 
FJBC. Let us note that the remaining constraints are imposing conditions on the reciprocal distances 
between two points (the centers of mass of the spheres) exactly as the constraints of the FHC given 
in Eq. ( !33|) do. As a consequence, it is possible to use the same strategy proposed in Subsection 13.11 



in order to fix the conditions of Eqs. 
noise distribution: 

N K 

1^pi{j^i,k}) 



Accordingly, we introduce the following non-gaussian 



-1 



=1 k= 



X exp 



1 

1 

AD 



T^^N,K exp 



N K-l 



=1 fc=i 



^Ikdt 



— / i^Nxdi 







(79) 
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The potential W{{Ri^k,Ui f^i^ t^2, f^s) is a function of the solutions {Ri^k,u^k} ^'i^- (IZ3^ZH]) and 
thus of the noises Ui^^ for i = 1, . . . , N — 1 and k = 1, . . . , K — 1 and Vn,k- The form of W can 
be chosen following the same criteria [H - El discussed in Subsection 13.11 A slight difference from 
the FHC is that, due to the splitting of the constraints (1^ into the three sets of equations given 
in ( 172]47il) . the FJBC potential W{{Ri^k,Ui ,,}, Pi, (^2, /^s) will depend on three parameters /32 and 
/^s instead of one. The underlying idea is however the same. When these parameters approach 
infinity, W{{Ri^k,un,}, (^2, /^s) should be chosen in such a way that it becomes infinite outside the 
region in the coordinate space in which the constraints ( 172147^ are satisfied and zero otherwise. Let's 
now derive W{{Ri^k,Ui f^i, (^2, f^s) explicitly. We have to require that the solutions Ri^k,Uik of 
Langevin equations (177]) and (!78l) satisfy the constraints (172147^1) . We recall that the Ri^k,uik^ are 
purely functions of the noises z/j fc, so that the relations (172147^1) are conditions on the noise degrees 
of freedom. Moreover, Eqs. (172H7^ are of the form: 

Ra — Rb = dab (80) 

In other words, as previously stressed, they constrain the distance between two points a i b exactly 
as the constraints of the FHC of Eq. (1351) . In order to implement constraints of this kind we can 
imagine, as it has been done in the previous Subsection, that the two points a and b are connected 
together by a spring with rest length dab > 0. In the limit of infinite elastic constant, the spring is 
frozen in its rest position and the distance between the points is fixed to dab- These consideration 
suggest that constraints of the kind given in Eq. flHOj) may be fixed using the elastic potential 

Vab{\Ra - Rb\) = ^ {\Ra " Rb\ ' dabf (81) 

There is however no reason for restricting ourselves to elastic interactions. Any two-body interaction 
that freezes the distance between two points in the limit of infinite strength is suitable. Of course, 
the related potential should have an infinitely sharp peak in the region for which \Ra — Rb\ = dab- 
Moreover, it has to be of the form V{\Ra — Rb\) = P{\Ra — Rb\)i where / is an arbitrary function, 
otherwise it will no longer be possible to use the identity fl^U]) . For example, a valid potential is the 
following: 

V:,,\Ra -R,\) = P{\Ra- Rb? - dlY (82) 

It is easy to show that the above potential is completely equivalent to the elastic one given in Eq. flHT]) . 
This becomes clear if we rewrite V^'^(|J2a — Ri,\) in the form 

V:,i\Ra -Rt\)=Pi\Ra-Rb\- dabf i\Ra - Rb\ + dabf (83) 

Of course, the minimum of V^f^{\Ra — Rb\ at \Ra — Rb\ + dab = will never be reached because 
dab > 0. Therefore, exactly as in the case of the elastic potential ( 15T|) . V^jj{\Ra — Rb\) imposes only 
the condition \Ra — Rb\ — dab = 0. Moreover, it is possible to check that V^'^(|i?a — Rb\) verifies all 
conditions 1. - 3. of Subsection 13.11 In particular, in the limit /3 = +oo V^'^d-Ra — -R^l) is zero if the 
constraints are satisfied and infinite in the opposite case. The advantage of choosing the potential 
of Eq. (!82|) is that it leads directly to a delta function in the desired form 6{\Ra — Rbl"^ — d^i^), so 
that we do not have to pass through all the intermediate steps performed in Eqs. (1621) and (1631) in 
order to arrive at the final generating functional ZpHc[J] of the FHC given in Eq. (I64p . Taking into 
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account all the above considerations, to implement the constraints of Eqs. (172147^ we will choose the 
following potentials: 



Vi,k{\ri+i,i - Vi^kl) = A//32 



\K-k\ 



I 
k 



N -I 
K-1 



15) 



Vk{\rN,K - rN,k\) = A//33 



rN,K - rN,k 


2 


K-k 




a? 









k 



K-1 



Here (3i, (32, /Ss are real parameters and are supposed to be very large. Eventually, we will take the 
limit Pi, (32, (3^ oo. Let us note the change of length scale with respect to the FHC Here the 
smallest scale is the distance Al between the centers of mass of two contiguous spheres belonging to 
the same bar. As a consequence, the analog of Eq. (fT2!) in the present case is: 



K 



K 



VpjBcm = J2AlVE{R^,k{t)) + Yl ^l^Vi{R,,k{t),R,,k'{t)) 



17) 



k=l 



k,k'=l 



The above form of the potential guarantees the correct passage to the limit AZ — )■ and K — +oo 
in which the bars become continuous and one- dimensional systems with uniform mass distribution. 
Indeed, it may be easily checked that in this continuous limit Eq. (187|) reduces to the potential for a 
chain of one- dimensional bars with uniform and continuous mass distribution given in Eq. fl28p . 

Going back to the main problem of constructing the generating functional of the FJBC, we define 
W{{Ri^k}) as the following linear combination of the potentials 
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(88) 



By substituting the right hand side of Eq. (188|) in Eq. (1791) . we obtain the noise distribution 
'Dp{{i'i k})- Knowing the noise distribution it is possible to construct the generating functional of 
the FJBC in path integral form: 
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At this point, in analogy with the steps f H71 - H9|) made in Section [STTI we insert in the generating 
functional ( 15^ the quantity: 



N K-1 



i=i k=i 



T^'l'N,KS{fN,K — Rn,K,vn,i 



(90) 



Clearly J = 1, so its insertion in Eq. ( 189|) does not change the physics of the problem. It is also 
possible to check that: 



N A'-l 



nndet 



d_ 
di 



Vri^kSivi^k - V. 



i,k 



det 



d_ 
di 



T>rN,K5{rN,K - 'i^n,k) (91) 



_i=i k=i 

As a result of the insertion of I, we get after an easy integration over the noises Ui^^ and iyN,K- 
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(92) 
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Finally, we apply Eq. (160|) in order to perform the limits P2, ^ +00. The result is: 
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(93) 



This is the expression of the generating functional of the discrete FJBC in the shish kebab approxi- 
mation. 



4 The double continuous limit of the FJBC 

Before performing the continuous limit of the generating functional (l93ll . it is important to recall 
how the FJBC is constructed. The FJBC consists in one-dimensional systems with an uniform 
and continuous distribution of mass as it was explained when discussing the classical case and in 
Subsection 13.21 In our approach each bar has been approximated by a shish-kebab model, in which 
K small spheres are put together in a row by imposing the constraints f l72] - [71|) . The distance between 
the centers of mass of the spheres is given by the quantity A/ of Eq. fl^TI) . One can see the spheres as 
small masses m localized at regular intervals on a segment of length a. The three dimensionality of 
the spheres is adjusted by tuning the diffusion constant D. In the present case D has been choosen in 
such a way that it coincides with the diffusion constant of a small sphere of diameter A/. Of course, 
if the length a of the bar is finite, the more spheres we will use, the better will be the approximation 
of the bar. To reproduce the FJBC, an infinite number of spheres is necessary. Indeed, when K is 
infinite, the diameter of the spheres approaches zero as it is shown by Eq. fl67|) and the shish-kebab 
models coincides with an uniform distribution of point-like masses distributed on a segment of length 
a. This is exactly the one-dimensional bar that has been taken as the basic unit in the FJBC model. 

As a warming up exercise, the probability function of the FJBC will be computed in the limit 
A/ — and K — )■ +00. As explained before, this limit corresponds to the one-dimensional chain with 
uniform mass distribution without the shish-kebab approximation. In that limit L and L coincide 
and the length of each bar will be equal to a because the decrease of the length due to the way in 
which the bars are joined together becomes negligible when A/ — )■ and K -\-oo. 

First of all, we derive the expression of the mass m of a single sphere in terms of two macroscopic 
parameters, namely the total length of the chain L and its total mass M. To compute M, we subtract 
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from the total mass NKm of bars the mass {N — l)m of the A^ — 1 spheres which disappear because 
they are identified at the junctions between the bars. The resuh is M = NKm — {N — l)m, Thus 
we may write: 

M , , 

m = (94) 



Using Eq. ( I7T1) . we get: 



N{K - 1) + 1 

Ma M 
-- = ^Al 



KL L 



where is the uniform mass density. The analog of Eq. ( l66ll in the FJBC case is: 



1 



1 ^ A, 



AD 2kTT 2L 

Next, we provide the prescriptions to pass to the continuous bar limit A/ — )■ and K — )■ +oo: 

K-l ,n K-l 



(95) 



(96) 



k=l "^0 fc,fc'=l "^0 



dldl' 



ri,k{t) ri(t,l) 
r,(t,0) = r,,i(t) 



ri{t,a) = r"i+i,i(t) 



Using the above formulas (!96] - ll00p is is not difficult to show that: 



7K=oo 
'FJBC 



[J] 



lim 

/3i,/32,/33->+oo 



where 



X exp 
X exp 



n / '^^iitj) 
1=1 ^ 

.•1 Jo 



i=l 



exp 



exp 



-Y.!'dtif\t) 

i=i 

N-l i-t; 



/33 / dtlf\t) 







dl 



1 M 



r^{t,l) + Ji{t,l)-ri{t,l) 



dldl 



dl 



2kTT 2L 

, /|ri(t,/)-n(t,/')P 



I/-/' 



|n+i(t,0)-n(t,/)p \a-l 



\2\ 2 



|2\ 2 



dl 



|rAr(t, a) — rAr(t, /)p \a — l 



2\ 2 



(97) 

(98) 
(99) 

(100) 



(101) 

(102) 
(103) 
(104) 
(105) 
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By taking the limit /32, /^s — )■ +00 we finally obtain the probability function of the FJBC: 



zK=oo 
'FJBC 



[J] 



N 

n J ^^^(^'0 



i=l 
■ TV 



exp 



N 

E 

i=l 



i=l 



\ri{t,l) - ri{t,l 



/\12 



if 

dt 
I -I 



dl 



M 



2kTT2L 

N-l 
i=l 



r^{t,l) + Ji{t,l)-ri{t,l) 



la — l\ 



X 6 



\rN{t,a) — rN{t,l)\'^ \a — l\ 



(106) 



We would like to check if, in the limit of a continuous chain (E]), the generating functional of the 
FJBC reduces to the partition function of the GNLaM of Eq. fj65|) as it should be expected. Our 
starting point will be the generating functional Zpj'^lJ] in the version of Eq. fllOll) . First of all we 
expand the functions lf\t), . . . , I^\t) given in Eqs. (I102H105]) at the leading order in a. Since at the 
end the limit a — )• will be taken, there is no need to compute higher order terms. A straightforward 
application of the expansion in Taylor series gives for ll'^\t): 



M 



V 2kTT 2L 



rtl{t) + J^At)■ri,l{t)]a 



(107) 



To write down the above equation we have used the prescription (l99ll according to which rj(t,0) = 
rj^i(t). Let us consider now the term lf'\t)x. Due to the fact that we are working in the region of 
small values of a, it is possible to make the following approximation: 



ri(t,/) -n(t,/') 



ri{t,a) - ri{t,a) 



(108) 



drijtfl) 



dl 



Both members of the above equation coincide in fact with the modulus of the derivative 
up to higher order terms in the infinitesimal quantities I, I' and a. Using Eq. f llOSp to simplify the 
expression of Ii^\t) given in f llOSp . we obtain: 



\ri+i^i(t)-ri^l{t)\"^ 

|rjv(f,a)-rjv(t,0)| = 
a2 



1 A for i = 1 A^- 1 



for i = N 



-1 A 



(109) 



where A = dldl' ^^ V . A direct calculation shows that: 



A 



15 



(110) 



Similar calculations allow the computation of Ii'^\t) and l'^'{t): 



|r-,+i,i(t) -r,,i(t)p 



a 
5 



A^ 



|rAr(t,a) - rAr(t,0)p 



2 

11 ^ 



(111) 

(112) 
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As we see from Eqs. fll07H112"|) . in the action of the generating functional Zpj'^lJ] many degrees of 
freedom disappear when a becomes smalL The only remaining degrees of freedom are the variables 
— ''"i(^)O) for i = 1,...,A^ and rj^(t,a). All the other degrees of freedom ri(t,l) for / 7^ 
and I 7^ a can be simply integrated out in Eq. fllOip . because they do not appear in the functions 
Ii^\t), . . . , I]^\t) of Eqs. ( I102H105]) and thus in the generating functional Zpj'^lJ]. Their integra- 
tion will result in a constant overall factor C multiplied with the rest of the generating functional. 
Moreover, when a goes to zero, all terms containing positive powers of a vanish identically if these 
powers are not absorbed by an appropriate number of sums over all bars according to the prescription 
J2iLi ^ ^ lo particular it is easy to check that 



lim > 



i=l 



dtl 



(1) 



and 



lim 

a-i-O 



if 



dtl^^^ = 



(113) 



(114) 



As a consequence, for small values of a the leading order contribution to the generating functional 
^fTbcI-J] is provided by: 
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lim C 

2^+00 
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n / ^^^-iw 


exp 


j=i 





N j. 

i=i i 



X exp 



i=l 



dt 



- 1 



:ii5) 



Now it is possible to proceed with the continuous limit ([6]) following the prescription of Subsection 12.11 
The result is: 

7K=oo r 71 _ v'r 71 ''116) 



lim Z'^y^ciA = ZV] 

a^O,N^+oo,Na=L 



where 



Z'lJ] 



lim C / Vr^t, s) exp 



r dt [ ds {^§:r\t, s) + J{t, s) ■ r{t, s)) 
Jo Jo 

2 



X exp 



' 5 



dt 



ds 



dr{t, s) 



ds 



;ii7) 



Performing also the limit for infinitely large values of f32, we obtain up to an irrelevant proportionality 
constant: 



zy] 



X 6 



Vrit, s) exp 



dr{t, s) 



dt 



ds{^lir\t,s) + J{t,s)-r{t,s)) 



ds 



:ii8) 



This is exactly the partition function of the GNLcrM of Eq. (1651) as expected 
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5 Conclusions 



In this work an approach to the dynamics of an inextensible chain has been proposed. At the basis of 
that approach there is the observation that in a constrained mechanical system fluctuating in some 
viscous medium at constant temperature, the random noise acting on the system is subjected to 
constraints too. Following this simple observation, the constraints have been imposed here directly 
on the noise degrees of freedom. In this way the complications of having to deal with an extended 
set of degrees of freedom (the Lagrange multipliers) or with generalized coordinates, are absent. In 
the generating functional of the FHC introduced in Eq. fH6|) . the constraints have been fixed by 
means of potentials. We recall at this point that even in the case in which the strength of the 
potentials imposing the inextensibility constraints are very large, numerical simulations show that 
the fast modes related to the changes of the lengths of the basic units still dominate over the slow 
modes associated to the conformational changes of the chain [51 EI]. To eliminate this problem, in 
the generating functional of the FHC the limit in which the strengths of the potential (jH]) becomes 
infinite has been performed. The final expression of the generating functional of the FHC in this 
limit is shown in Eq. ( l64ll . 

In this work it has also been explored the strategy of Ref. [5], in which it has been proposed 
that an inextensible chain may be realized gluing together a set of inextensible basic units like 
for instance one-dimensional bars with uniform mass distribution. In Subsection 13.21 this strategy 
has been explored for a chain in which the bars have been replaced by the shish-kebabs displayed 
in Fig. [21 Each shish-kebab consists of K spheres that are aligned together in a row by means 
of suitable constraints in order to form the approximate shape of a bar. The final expression of 
the generating functional ZpjBc[J] for this type of chain may be found in Eq. fl^^ . In the limit 
in which the distance A/ between the centers of mass of the spheres vanishes identically and the 
number K of spheres becomes infinitely large while the product K/S.I remains finite, the shish-kebabs 
become one-dimensional bars and one obtains the model of a chain which has been called here 
FJBC. Its generating functional in path integral form has been given in Eq. fll06p . As a check of 
the consistency of our approach it has been verified that in the limit of a continuous chain both 
generating functionals of the FHC and FJBC coincide. Indeed, after performing the continuous 
limit, one obtains respectively the generating functionals of Eqs. fl65|) and flllSp that describe the 
same theory, namely the GNLaM of Ref. [29]. The consistency of the approach presented in this 
work in order to tackle the problem of the dynamics of an inextensible chain is also confirmed by 
the results of Ref. [32] • In has been shown there that the generating functional of the FHC given in 
Eq. ( 16^ consists of a statistical sum over all particle trajectories that satisfy the Langevin equations 
( 132|) together with the inextensibility constraints ( 133|1 as it is expected. In [30] it has been verified 
that in the continuous case the statistical sum in the GNLuM contains only the chain conformations 
for which R(t,s) = vitjs) and \dR{t,s) /ds\ = 1. These equations are respectively the continuous 
counterparts of Eqs. (13 2 p and (I33p . Additionally, we have explored in this work the Fokker-Planck 
formulation of the dynamics of the FHC in Eqs. (|55]) - (15^ . 

Finally, let us mention some issues that are still open. First of all, the hydrodynamic and self- 
avoiding forces are missing in our approach. To describe the dynamics of realistic chains these 
interactions should be taken into account. For this reason, up to now the applications of the GNLaM 
are limited to situations in which the speed of the elements of the chain is low, so that the hydro- 
dynamics interactions can be neglected. Let us note that the addition of nonlinear self-interactions 
acting on the beads of a set of coupled oscillators like that described by Eq. ( 150|1 leads to interesting 
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phenomena as described in Refs. [391 SHI IH] , where such kind of systems has been studied. Another 
important issue is to find suitable approximations that may be used to evaluate field path integrals in 
the presence of functional Dirac delta functions and nontrivial boundary conditions. This is exactly 
the case of the generating functional of Eq. fl65l) . So far, the partition function and the so-called 
dynamic structure factor have been derived in the semiclassical approximation [29| [30] . One could 
simplify the theory even further if it would be possible to extend the well known approximation of 
the Dirac delta function with a gaussian function also to the functional Dirac delta function imposing 
the inextensible constraints in the GNLcrM. At the moment, this extension has been achieved only in 
the case of the statistical mechanics of the FHC [26]. The generalization of this result to dynamics 
is still work in progress. 
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Appendix A — Proof of Eq. (1601) 

We start from the quantity 



A{i3) = C{(3) exp 



a/3 / dtfir,it),r,it)) 



z,j = l,...,N 



(119) 



where C{(3) is a normalization constant. We discretize the interval of time [0,t/] into M smaller 
intervals of length b. Of course, Mb = tf. is recovered in the continuous limit b — > and 

M — > +00: 



A(p) = lim Cm exp 



M 



-a(3 ^ bf{ri{tm), Vjitr 



Here = mb is the discretized time variable and C(/3) = lim Let us put 



(120) 



and 



Thus 



- = aPb 
e 



M 



c^^(^) = n ttV 



A{P)= lim 

M— )- + oo 



TT — — ^exp 



f{ri{t„^),rj{tm)) 



We consider at this point the limit (3 — )■ +00 of A{f3): 

M 



lim A(p) = lim lim 



f (^i (^m) ; ''^j ('^m) ) 



'1211 



(122) 



(123) 



(124) 
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Permuting the limit e — with the hmit M — )■ +00 and 6 — 0, we get: 



hm AiB) = hm 

i\/-> + CX) 



M 



TT hm r exp 



m=l 



(125) 



When /3 starts to be very large and simultaneously e becomes very small, it is possible to use the 
Gaussian representation of the Dirac delta function 6(x) = lim^^Q rC ^ and to write: 

(47re)2 



lim A(p) = lim 

J\-/— > + oo 



M 



J]5(/(n(t„),r,(U)) 



m=l 



(126) 



Here we have exploited the fact that, if the limit for e — of a function (^^(e) exists for every 
m = 1, . . . ,M, i. e. lim,^,^_+o 5'm(e) = gm{0) with g^iO) ^ and \gj^)\ < +00, then we may write 
the following relations: 



5'i(0)5'2(0) ■ ■ ■ guiO) = lim fi'i(ei) lim 5(2(^2) ■ ■ ■ hm gui^^M] 

= lim(^i(e)^2(e) ■ ■ ■^M(e)) 

e->-0 



(127) 



The right hand side of Eq. f ll26p is the definition of the functional delta function 6{f{ri(t),rj(t))) 
which, at each instant t, is concentrated in the points for which f{ri{t), rj{t)) = 0. In other words: 



lim A{(3) = 6{f{r,{t),r,m 



(128) 



Remembering the definition of A{I3) of Eq. (I119p . we obtain: 



lim C (/3) exp 

/3— s>+oo 



af3 / dtf{r,{t),r,{t)) 



5{f{n{t),rM) 



(129) 



This proves Eq. (160|) apart from the infinite normalization constant C(oo). 

Alternatively, one can show Eq. fl60l) using the Fourier representation of the Dirac delta function: 



(130) 



5{f{r,{t,s),rAt,s)))= lim / pA(t, s)e-*/o^ '^t/o^ '^«A{t,s)/(r..(M),r,(M))e-/J^ /o^ 
Performing the simple gaussian integral in A(t, s) we obtain: 

5(/(n(t,s),r,(t,s))) = lim /"l?A'(t,s)e-^o^'^*^o^'^«^^e-^^o''^*^o^'^-^'("'(*'^)'''^(*'^)) (131) 



where A' is the shifted field A'(t,s) = A(it, s) + iaf{ri{t, s),rj{t, s)). Thus, apart from an irrelevant 
constaint, we may write the relation: 



5(/(ri(t,s),r,-(t,s))) = lim g't io'^ '^*/o^'^^/'(--(*.^).-.(M)) 

a— >+oo 

which coincides exactly with Eq. (|60ll . 



(132) 
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